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Abstract: We consider chiral perturbation theory for SU(2) at finite temperature T in a constant 
magnetic background B. We compute the thermal mass of the pions and the pion decay constant to 
leading order in chiral perturbation theory in the presence of the magnetic field. The magnetic field 
gives rise to a splitting between M n o and M n ± as well as between F n o and F n ± . We also calculate 
the free energy and the quark condensate to next-to-leading order in chiral perturbation theory. 
Both the pion decay constants and the quark condensate are decreasing slower as a function of 
temperature as compared to the case with vanishing magnetic field. The latter result suggests that 
the critical temperature T c for the chiral transition is larger in the presence of a constant magnetic 
field. The increase of T c as a function of B is in agreement with most model calculations but in 
disagreement with recent lattice calculations. 
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1 Introduction 

While the QCD Lagrangian in the limit of zero quark masses has chiral symmetry, the true ground 
state of QCD does not respect this symmetry. Chiral symmetry is broken spontaneously in the 
vacuum by quantum effects. Specifically, the SU(Nf)i, x SU(Nf)n symmetry of £qcd is broken 
down to SU(Nf)v, and according to Goldstone's theorem this gives rise to as many massless spin-zero 
particles as there are broken generators. For Nf fermions, this implies Nj — 1 massless excitations 
and in phenomenological, we have applications Nf = 2 or Nf = 3. The symmetry breaking is 
apparent in the low-energy spectrum of QCD from the existence of the three very light pions. The 
fact that the pions are not strictly massless is due to the fact that chiral symmetry of nature is only 
approximate. The associated explicit symmetry breaking in QCD is built in by finite quark masses 
in £qcd- 

Lattice simulations are currently the only way to calculate the properties of QCD from first 
principles. However, due to the sign problem, one is restricted to small values of the baryon chemical 
potential. Other approaches involve low-energy effective theories or model calculations that share 
some of the properties of QCD such as the same pattern of symmetry breaking. Examples are 
chiral perturbation theory (ChPT), the Nambu-Jona-Lasinio (NJL) model, linear sigma models, 
and quark-meson (QM) models. At long wavelengths where the relevant degrees of freedom are the 



almost massless Goldstone bosons, one would naturally like to have an effective low-energy theory 
that contains only these degrees of freedom. Such an effective field theory may simplify calculations 
of the long-distance properties of the system and one can obtain model-independent predictions with 
a minimum of assumptions. 

An example of a low-energy effective theory for QCD is chiral perturbation theory [1—4]. The 
chiral Lagrangian that describes the (pseudo) Goldstone bosons is uniquely determined by the global 
symmetries of QCD and the assumption of chiral symmetry breaking. The Lagrangian C c s consists 
of a string of operators that involve an increasing number of derivatives or quark mass factors, each 
multiplied by a low-energy constant (LEC) l{. However, QCD is a confining and strongly interacting 
theory at low energies. Thus the coefficients l{ of the chiral Lagrangian cannot be calculated directly 
from QCD. Instead, the coefficients of the effective Lagrangian are fixed by experiments. 

Chiral perturbation theory provides a systematic framework at low energies. It is not an expansion 
in powers of some small coupling constant, but it is a a systematic expansion in powers of momenta 
p where a derivative counts as one power and the quark masses count as two powers [2]. Chiral 
perturbation theory is a nonrenormalizable quantum field theory in the old sense of the word. This 
means that a calculation at a given order n in momentum p, requires that one adds higher-order 
operators in order to cancel the divergences that arise in the calculations at order n. This implies 
that one needs more and more couplings li as one goes to higher loop orders, and therefore more 
experiments to determine them. However, this poses no problem, as long as one is content with 
finite precision. This is is the essence of effective field theory [5-7]. The order p 4 -divergences in 
ChPT were determined in Ref. [2], while the order-p 6 divergences were calculated in Ref. [4]. 

The thermodynamics of a pion gas using ChPT was studied in detail in a series of papers 25 years 
ago [8-10]. The temperature-dependence of the pion mass M n and the pion decay constant F n were 
calculated to leading order (LO), while the the pressure and the temperature dependence of the 
quark condensate were calculated to next-to-next-to-leading order (NNLO) in chiral perturbation 
theory. In the chiral limit, there are two scales at finite temperature, namely the pion decay constant 
and T, and chiral perturbation theory is then an expansion in powers of T 2 /F%. In the present paper, 
we give the calculational details of a number of quantities in the presence of a constant magnetic field 
B. These include M n o, M n ±, F w o, F n ± at leading order and the free energy and the (normalized) 
quark condensate at next-to- leading order (NLO). 

QCD in external magnetic fields has received a lot of attention in recent years. This is not a purely 
academic question since the properties of QCD in strong magnetic fields B is relevant in several 
situations. For example, large magnetic, B ~ 10 14 — 10 15 Gauss, exist inside magnetars [11, 12]. If 
the density of the core is sufficiently high, it contains quark matter. In that case, the core may be 
color superconducting and so it is important to study the effects of external magnetic fields in this 
phase [13-20] . Similarly, it has been suggested that strong magnetic fields are created in (noncentral) 
heavy-ion collisions at the Relativistic Heavy-Ion Collider (RHIC) and the Large Hadron Collider 
(LHC) and that these play an important role [21-23]. In this case, the magnetic field strength has 
been estimated to be up to B ~ 10 19 Gauss, which corresponds to \qB\ ~ 6m^, where \q\ is the charge 
of the pion. Even larger fields could be reached due to the effects of event-by-event fluctuations, see 
for example [24]. 

This has spurred the interest in studying QCD in external fields. At zero baryon chemical potential 



this can be done from first principles using lattice simulations. Recently, lattice simulations in a 
constant magnetic background B have been carried out [25-27]. At finite \ib this is very difficult 
due to the infamous sign problem. Therefore one often resorts to effective theories that share some 
of the features of QCD, such as chiral symmetry breaking. 

A low-energy effective theory that provides a systematic framework for systematic calculations is 
chiral perturbation theory. Chiral perturbation theory has been used to study the quark condensate 
in strong magnetic fields at zero temperature [28-31] and finite temperature [32]. In Ref. [33], the 
thermal corrections to M n o and F w o were computed, while the quark-hadron phase transition was 
investigated in Ref. [34]. The effects of external magnetic fields on the chiral transition have been 
studied in detail using the NJL model [35-44], the Polyakov-loop extended NJL model [45, 46], the 
QM model [43, 44, 47-49], the linear sigma model [50], the (P)QM model [51, 52], and the MIT bag 
model [53]. 

In the present paper, we study pions at finite temperature in a constant external magnetic field B 
using chiral perturbation theory and it is organized as follows. In Sec. II, we briefly discuss ChPT 
in an external magnetic field. In Sec. Ill, we calculate the leading correction to the pion mass as 
well as the pion decay constant. In Sec. IV, we calculate the free energy and the quark condensate 
to next-to-leading order in chiral perturbation theory. In Sec. V the numerical results are presented 
and we summarize in Sec. VI. The necessary sum-integrals are listed in Appendix A, while explicit 
calculations of those sum-integrals are presented in Appendix B. 

2 Chiral perturbation theory 

As explained in the introduction, chiral perturbation theory is a low-energy effective field theory 
that can be used to systematically calculate physical quantities as a power series in momentum p. 
The effective Lagrangian is given by an infinite string of operators involving an increasing number 
of derivatives or quark masses. Schematically, we can write 



where the superscript indicates the powers of momentum. In Euclidean space, the leading term is 
given by 



The first term is the Lagrangian of the nonlinear sigma model and the second term is the kinetic 
term for the gt Here U = exp[ir a 7r a /-F] is a unitary SU(2) matrix, where n a are the pion fields and 
T a are the Pauli spin matrices. The low-energy constants M and F are the tree-level values for the 
pion mass M n and the pion decay constant F n , respectively. Moreover is the covariant derivative 
which replaces the partial derivative: 
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where Vu is an external vector field and is an external axial vector field. The next-to-leading 
order terms in the low-energy expansion are [2, 10, 31] 1 
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Here F 1 ^ = d^F* - d u F 1 ^ - Fj], where l = L,R, and F* = v^ + a^ and F% =v )JL -a fl . 
The Lagrangian 

£(6) 

is very complicated as it contains more than 50 terms for SU(2) [4]. However, 
only one term is relevant for the calculations [29-31], namely 
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We consider in the following a constant external magnetic field so is an Abelian gauge field, 
= ^qA^T 3 , and = 0. Here q is the electric charge of the pion and we choose the four-vector 
potential to be A^ = B8a2 x X- By expanding the Lagrangian Eq. (2.2) to fourth order in the pion 
field 7r, we obtain 

£(2) = _ F 2 M 2 + 1 ^0)2 + 1 M 2 ^0)2 + ( ^ + iqA ^+ ^ _ iqA ^~ + M 2 7T + 7T~ 

-2 [(vr ) 2 + JJtt+tt] (a M 7r+)(a M 7r-) + (vt + tt" )] 2 1 , (2.6) 

where we have defined the complex pion fields as n± = (tti ± m^j In the same manner, we can 
expand C^: 

£(4) = 7*% + !§(^) 2 7r+7T- + ^qFp, [(d^-)(d u -K+) + MjAOr+O] + (/a + i*)^ (vr ) 2 



F 2 



F 2 



M 4 



+2(l 3 + / 4 )9 2 - 7r+7r + Z 4 : i 2 -( a M 7r0 ) 2 + 2Z 4 : iT(5 M + ^)7r + (^-^)7r . 



F 2 



F 2 



F 2 



(2.7) 



We note in passing that our expression (2.6) for the truncated Lagrangian £^ differs from the 
expression found in Refs. [29, 30] since they use a different parametrization for the unitary matrix 
U, namely the Weinberg parametrization. However, we obtain the same expressions for physical 
quantities independent of parametrization [54]. This is simply the statement that physical quantities 
are independent of the coordinate system used. 



1 The last two operators of Eq. (2.4) are so-called contact terms. The term —hiM 4 contributes to the vacuum 
energy. Loop corrections to the vacuum energy are divergent and renormalizing the couplings such as hi eliminates 
(some of) the divergences and is needed to render the vacuum energy finite. We will not need them in the following 
since we are subtracting the divergences in the vacuum diagrams that appear for B = 0. The remaining divergences 
are then eliminated by renormalizing the low-energy couplings U. 



At this point it is appropriate to briefly discuss the number of Goldstone modes in the presence of 
an external electromagnetic field. Due to the different electric charges of the u and d quarks, flavor 
symmetry is broken in an external electromagnetic field. In particular, the axial S77(2)^ -symmetry 
is broken down to a [/(l)^-symmetry, which corresponds to a rotation of the u and d-quarks with 
opposite angles. The formation of a quark condensate breaks this Abelian symmetry which gives 
rise to a Goldstone boson - the neutral pion ir° [28] . The charged pions are thus no longer Goldstone 
modes. In fact, the presence of the external field allows for an effective mass term even in the chiral 
limit when M = 0, namely the first two terms in Eq. (2.7). 

The chiral Lagrangian comes with a number of undetermined parameters or low-energy constants 
(LECs). These parameters can be determined by experiments. However, loop corrections involve 
renormalization of these parameters and the physical quantities are no longer equal to the parameters 
of the chiral Lagrangian. For example, F and M are no longer the measured pion decay constant and 
the physical pion mass, respectively. The relation between the bare and renormalized parameters 
has been found in [2] and can be expressed as 2 
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where 7$ are coefficients that are tabulated in [2] and li are scale-independent parameters, i. e. they 
are the renormalized running couplings evaluated at the renormalization scale A = M. 

In the present calculations, we need 73 = — |, 74 = 2, 75 = — |, and 76 = — | [2, 3]. The coupling 
c 34 m 

£(6) 

can be renormalized as in [4] by using 
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We close this section by briefly discussion the bare propagators of the pions in a magnetic backround. 
The classical solutions to the Klein-Gordon equation in a constant magnetic field B are well known 
and the dispersion relation is given by 

(Et, p f =p 2 z +M 2 + (2m + l)\qB\ , (2.10) 

where m = 0, 1, 2, ... denotes the mth Landau level, q is the electric charge of the pion, and p z is the 
spatial momentum in the z-direction. The subscript ± denotes tt^ and we note that the dispersion 
relations for the two charged pions are identical. The imaginary-time propagator for a charged pions 
is then A 7r ±(Po ) p) = VC^o + Pz + M 2 + (2m + l)|gP|), while for neutral pions it is given by the 
usual A 7r o(Po J p) = l/(^o + P 2 + A^ 2 ); where Pq = 2nirT is the nth Matsubara frequency. 



3 Pion masses and pion decay constants 

In this section, we discuss the thermal pion masses and the pion decay constants. In order to 
calculate the thermal pion masses, we consider the inverse propagators at one- loop. For example, 
the inverse propagator for the neutral pion can be written as 

A^PcP) = P 2 + M 2 +S 7ro (P ,p) , (3.1) 

2 In contrast to ref. [2], we have introduced the renormalization scale A in the definition of the sum- integrals instead 
of as a part of the relation between the bare and renormalized parameters. See Appendix A for details. 



where the self-energy is given by 

S^(P , P) = [Ao - AZ ]P 2 + [B - AZ ]M 2 , 



(3.2) 



where AZq is the wavefunction renormalization counterterm for ir°. and ir^, respectively. The 
coefficients Aq and Bq are given by 
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The sum-integrals and ^Qf are defined in Appendix A. The counterterms AZq is chosen such that 
the self-energy Yl wo is independent of the momentum of the pions [55]. This implies that AZq = Aq. 
The propagator then reduces to 
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Inserting Eqs. (3.3) and (3.4) into Eq. (3.5), we can then read off the mass m^o from the propgator 
and we find 
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Renormalization of I5 and are carried out according to Eq. (2.8), which shows that we can replace 
2Z 5 - l 6 by (l 6 - l 5 )/6(4ir) 2 [29, 30]. After renormalization, one obtains 
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where the function Jb(M) is defined by 
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and x = 2 j^L . In particular, note that -fe(O) = |gi?|log2. The thermal integrals J n ((3M) and 
J^{j3M) are defined in Appendix A, and evaluated at e = 0. Moreover, C ( — &>2/) = ^C( x ~ 



b, y)\ x =o, where £(s,x) is the Hurwitz zeta-function and b is a real number. The physical pion mass 
M 2 in the vacuum is given by 
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The result (3.8) is in agreement with the calculations of Ref. [33]. In the limit B — > 0, Eqs. (3.8)- 
(3.9) reduce to the result of Ref. [8]. We notice that temperature dependence of the charged pion 
mass M w ± is the same as for vanishing magnetic field since the only loop correction to the mass 
involves a neutral pion, cf. Eq. (3.7). The only difference is a temperature- independent constant 
proportional to (qB) 2 . This constant survives the the chiral limit, M — > 0, which shows that the 
charged pion is no longer a Goldstone mode in the presence of an electromagnetic field. 

We next discuss the pion decay constant F w o for the neutral pion. The components of the axial 
current A^ are given by [54, 55] 
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In order to calculate the pion decay constant, we need to evaluate the matrix element F w o = 
(0|^4^|-7r ). The matrix elements are proportional to ip^. Moreover, in the leading term —Fd^TT 
of Eq. (3.12), we must take into account wavefunction renormalization of the pion field given by 
Eqs. (3.3). A straightforward calculation of the matrix element then gives 
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In a similar way we can calculate the pion decay constant for the charged pion, F n ± , and one finds 
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After renormalization, we find 
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where the thermal integrals J\((3M) and Ji(j3M) are evaluated at e 
constant F n in the vacuum is 
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The result (3.15) agrees with the calculations of Agasian and Shushpanov [33]. In the limit B — > 0, 
F n o = F n ± and they reduce to the result of [9] as they should. At zero magnetic field, the decay 
constants are identical F n o = F w ± but there are two pion decay constants at finite temperature; one 
for the time component and one for the spatial component of A^. The difference between them is 
an order-p effect [56] and this the calculation of the diffference is beyond the scope of this paper. 



4 Free energy and quark condensate 



In this section, we compute the free energy and the quark condensate to NLO in chiral perturbation 
theory. We are interested in the contributions to the free energy that are due to a nonzero magnetic 
field and finite temperature. We therefore write the contribution to the free energy at the nth loop 
order, J- n , as a sum of three different terms 

J~ n = T n + J- n + J- n , (4-1) 

where J~^ c is the contribution in the vacuum, i. e. for B = T = 0, Jjf is the zero-temperature 
contribution due to a finite magnetic field, and J-£ is the finite-temperature contribution. The last 
term also depends on B. The strategy is to isolate the term J-^ ac and subtract it from Eq. (4.1). 
This term contains ultraviolet divergences which are removed by renormalization of the low-energy 
constants of the chiral Lagrangian in the usual way and the renormalized expression for J\^ c repre- 
sents the contribution to the vacuum energy of the theory from the nth loop order. The term J-^ 
generally contains additional ultraviolet divergences as well and it is rendered finite by the same 
renormalization procedure. The term J 7 ^ is finite. 
The free energy at tree level is given by the first term in Eqs. (2.6) and (2.7): 

Jo = I# 2 - F 2 M 2 . (4.2) 

The first term in Eq. (4.2) is needed when we renormalize the one-loop contribution to the free energy 
in the next subsection. The second term is independent of B, so we neglect it in the following. 

4.1 One-loop free energy 

The one-loop contribution T\ to the free energy can be written as the sum of the contributions from 
the charged and neutral pions, and tree-level terms from 
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= [Po +Pl + M 2 B ] , (4.5) 

The tree-level terms are necessary in the renormalization of the vacuum energy, but according to our 
strategy we do not need to consider these since they are independent of the magnetic field. Using 
the expressions for the various sum-integrals in Appendix A, we obtain 
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The divergence that depends on the external magnetic field B is removed by wavefunction renor- 
malization, i. e. by the replacement 
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in the tree-level term \B 2 in Eq. (4.2). Subtracting the expressions for the diagrams for T = B = 0, 
we remove the remaining divergences and obtain a finite result. One finds 
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4.2 Two- loop free energy 
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Figure 1. Feynman diagrams contributing to the free energy at order p 4 



The order-p 4 diagrams that contribute to the free energy are shown in Fig. 1. Dotted lines denote 
a neutral pion and solid lines denote a charged pion. Wiggly line denote the background field \qB\. 
The two-loop contribution Ti to the free energy can be written as 
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where the expressions for the diagrams are 3 
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3 It was shown in Refs. [29, 31] that the diagram J^c vanishes at T = 0. It also vanishes at finite temperature. 



Using the expressions for the sum-integrals in Appendix A, we obtain 
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After having subtracted the expressions for the diagrams for T = B = 0, we are still left with some 
simple poles in e. These remaining divergences are now elimated by renormalizing the couplings Z3, 
I5, and Iq according to the prescription (2.8). Since the renormalized couplings l\ are evaluated at 
the scale A 2 = M 2 , the final result simplies significantly. After a lengthy calculations, we obtain 
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Note in particular that all the divergent terms involving the thermal integrals J\ and Jf cancel 
amongst themselves and so we can evaluate the remaining ones at e = 0. 

If we express the contributions J-f and J-f in terms of the physical pion masses M n o(0) and 
M„.±(0) using Eqs. (3.8) and (3.9) instead of the parameterM, the dependence on £3, Z5, and Iq 



cancels in the sums J-f\ 2 and J~\+2- This yields 
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where x-t = M 2 ± (0)/2\qB\. The result for the vacuum energy J 7 j^ 2 ^ s i n agreement with the 
calculation of Ref. [31] 4 . Since M n ± is nonzero in the chiral limit (M = 0) and it appears in the 
thermal integrals Jq in Eq. (4.27), the NLO correction to the free energy does not vanish in contrast 
to the case of zero magnetic field. 

We next discuss the quark condensate in the presence of the the magnetic field B. At finite 
temperature, the quark condensate is [10] 



(qq) = (0\qq\0) 
where c is a constant defined by 
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where m is the quark mass. In the chiral limit, we have c = 1 [10]. In this case, the quark condensate 
reduces to 
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where we have used the recursion relations for J n and . We notice that we in the limit B — > 
recover the result by Gerber and Leutwyler [10] 



(4.31) 



where we have used Ji(0) = 47r 2 /3. 

4 Note that the final result of Ref. [31] is expressed in terms of the pion mass in the vacuum and so the terms 
involving Z5 and le are not absorbed in the order p 4 term. 



5 High-temperature expansion 



In this section, we discuss the expansion of our results in the limit of weak fields, i.e. for \qB\/T 2 <C 1. 
For simplicity, we restrict our analysis to the chiral limit, M = 0. 

For values \qB\/T 2 < 1, the sum over Landau levels in J B (f3M) converges very slowly. It is 
advantageous to use the Euler-McLaurin summation formula [57] to evaluate J B in this regime [50] . 
The sum over integers m can be written as 
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where B n are the Bernoulli numbers and {x} means the fractional part of x. The last term in 
Eq. (5.1) is the remainder. We will illustrate the use of Eq. (5.1) in the case of Jf. Setting n = z, 
k = 0, and I = oo, we obtain 
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where we have defined E(z, M, k, T) = \J z 2 + M 2 /T 2 + (2k + l)\qB\/T 2 . Integrating over k in the 
first term of Eq. (5.2), we obtain 
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where a 2 = \qB\/T 2 . The Euler-McLaurin formula is next used to expand the integrals Jjf in a 
Taylor series around B = 0. This is particularly simple in the chiral limit where Jjf only depends 
on the ratio \qB\/T 2 , and the series is in powers of y/qB/T. Integrating the first and the last term 
in Eq. (5.3) by parts, we can write 
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where the argument of J n is /3a. The expansion of the integrals J n (x) can be found in e.g. [60] and 
the leading terms are J\ = 4ir 2 /3 — Ana, J2 = 2n/a, and J3 = 7r/a 3 . This yields 
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Inserting this expression into Eqs. (3.15), and (3.16), we obtain in the chiral limit 
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In the same manner, we can expand the temperature-dependent part of the quark condensate 
around qB = 0. A straightforward calculation yields 



(qq) = (0\qq\0) 
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6 Numerical results and discussion 



(5.8) 



We need to evaluate the integrals Jq and Jf . By expanding the Bose-Einstein distribution function 
in Jq and J/ 3 , and integrating over three- momenta, we write it as a double sum involving a modified 
Bessel function K[{x) of order I = 0: 
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For values \qB\/T 2 > 1, the sums converge fast. In the remainder of the paper, we use \qB\ 
5(140MeV) 2 and so we use the expressions (6.1)-(6.2) in the numerical work. 



Thermal pion masses 



1 1 1 1 






Neutral pion mass 




Charged pion mass 










1,1, 



Figure 2. Pion masses M^ and M^ ± normalized to their vacuum valued M% as a function of T. 



In Fig. 2, we show the pion masses M 2 and M 2 ± given by Eqs. (3.8) and (3.9) normalized to 
their vacuum value M 2 function of T. We use the experimental value F n = 93 MeV. Moreover, 
the scale-independent difference of Iq and £5 is Iq — I5 = 3 ± 0.3 [58, 59]. In the remainder we use 
central value of this difference, i e. Iq — I5 = 3. We notice that the neutral pion mass is M 2 very 
large already at T = due the strong magnetic field \qB\ = 5(140MeV) 2 . 

In Fig. 3, we show the pion decay constant F® given by Eq. (3.15) at T = normalized by its 
vacuum value F n = 93 MeV as a function of \qB\/F 2 in the chiral limit as well as at the physical 
point. At T = 0, the dependence on the magnetic field is governed by the function Ib(M) defined 



in Eq. (3.10). In the chiral limit, F w o is a linear function of \qB\ since -Tb(O) = \qB\ log 2. In both 
cases, we note that pion decays constant is an increasing function of \qB\. We do not show since 
the functional dependence at T = is the same, cf. Eq. (3.16). 
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Figure 3. Pion decay constant F^o at T = scaled by its vacuum value F n as a function of \qB\/F%. 



In Fig. 4, we show the normalized pion decay constant F n o in the chiral limit as a function of 
T for \qB\ = (dashed line) and \qB\ = 5(140MeV) 2 (solid line). In vanishing magnetic field, the 
temperature dependence is given by 



F n o =F n [l 



T 2 \ 
12F 2 ) ■ 



(6.3) 



The pion decay constant is larger in a magnetic background. This is partly due to Ib(M) which is 
an increasing function of B and partly due to the fact that J B \qB\ < J\T 2 for all temperatures. 

In Fig. 5, we show the contribution F B to free energy in the chiral limit at one and two loops 
normalized to F^ as a function of \qB\/F 2 . We have omitted the term gT^jp log ^j^JJf as * s customary 
in the literature. The uncertainty in the coupling d = is large and even the sign is undetermined [58]. 
Varying d would give rise to large bands for the free energy as well as the quark condensate [31]. 
Since the experimental value is consistent with zero, we will simply use d = in the remainder of 
the paper. The two-loop result (4.26) for the free energy then reduces to the one-loop result (4.8) 
with an effective mass for the charged pion given by M 2 ± = ^^ypz (^6 — h)- 

In Fig. 6, we show the contribution F B to the free energy at the physical point at one - and 
two loops normalized to F^ as a function of \qB\/F 2 . The correction is tiny and the two-loop 
result starts to deviate from the one-loop result first at very large values of the magnetic field. 
Again the two-loop result is obtained from the one-loop formula (4.8) using the effective mass 
(^6 — ^5) instead of M 2 . 



M 2 ± =M 2 + 

In Fig. 7, we show the thermal part F T of the free energy as a function of T for qB = 5(140MeV) 2 
in the chiral limit normalized to the free energy Fq = — 7r 2 T 4 /30 of an ideal massless pion gas. In 
the chiral limit, the two-loop result (4.27) reduces to the one-loop (4.9) evaluated with the M 2 ± = 
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Figure 4. Normalized pion decay constant F„o for and \qB\ = and \qB\ = 5(140MeV) 2 as a function of T. 
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Figure 5. The contribution T B to the one and two- loop free energy scaled by in the chiral limit. 



(qB) 2 (lQ — 1$) / 3(4tt) 2 F 2 instead of M n = 0. For low temperatures, only the neutral pion contributes 
to the free energy since the contributions from the charged pions are Boltzmann suppressed. From T 
around 50 MeV onwards, the normalized free energy increases sharply and approaches the normalized 
Stefan-Boltzmann result 1 as T — > oo. It has the value of approximately 0.97 for T = 313 MeV, i. 
e. for T 2 /\qB\ = 1. The two- loop result lies below the one-loop result since in this approximation 
the charged pions are not massless but has an effective mass due to the magnetic field, M 2 ± = 
(qB) 2 (k-k)/3(4ir) 2 F 2 . 

In Fig. 8, we show the LO and NLO results for the normalized quark condensate in the chiral 
limit including the T = contribution. For comparison, we also show the LO and NLO results for 
the quark condensate for B = given by Eq. (4.31). In Fig. 9, we show the LO and NLO results for 
the normalized quark condensate in the chiral limit excluding the zero-temperature contributions in 
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Figure 6. The contribution J- B to the one and two-loop free energy scaled by at the physical point. 
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Figure 7. One and two-loop thermal part F T of the free energy normalized to the free energy Tq — — 7r 2 T 4 /30 
of an ideal massless pion gas. 

order to disentangle the effects of the magnetic field at T = and the finite-temperature effects. For 
comparison, we again show the LO and NLO results for the quark condensate for B = 0. Figs. 8 
and 9 show that the one-loop and two-loop results for the condensate are very close to each other in 
the entire temperature range. The LO and NLO results for the condensate are closer in a magnetic 
background than for B = and this suggests that chiral perturbation theory converges at least as 
well for finite B. 

Figs. 8 and 9 show that the quark condensate for nonzero B goes to zero slower than for vanishing 
magnetic field and the reason is two-fold. Firstly, there is an enhanced quark condensate at T = 0, 
which to leading order in ChPT is given by the function Ib(M). This is the familiar enhancement of 
the chiral condensate in a magnetic background. Secondly, there are finite-temperature effects. The 
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Figure 8. Temperature dependence of the quark condensate including the T = contribution normalized 
to its vacuum value \qB\ = 5 (140 MeV) 2 at LO and NLO in chiral perturbation theory. For comparison, we 
show the LO and NLO results for \qB\ = as well. 
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Figure 9. Temperature dependence of the quark condensate excluding the T — contribution normalized 
to its vacuum value \qB\ = 5 (140 MeV) 2 at LO and NLO in chiral perturbation theory. For comparison, we 
show the LO and NLO results for \qB\ = as well. 



basic mechanism is that J-j 8 is a decreasing function of the mag netic field B and so J{T 2 > Jf\qB\ 
for all nonzero B. 

Comparing the results for the condensate at B = and \qB\ = 5(140MeV) 2 , we see that its effects 
are quantitavely large. This is a due to a very strong magnetic field. For weaker magnetic fields, the 
difference between the two sets of curves will be smaller too. The results suggest that the critical 
temperature T c for the chiral transition is higher in a magnetic background. However, since ChPT 
is known to break down for large temperatures we do not make any quantitative statements about 



T c as a function of the magnetic field B. In Ref. [34], the authors are using chiral perturbation 
theory to investigate the quark-hadron phase transition as a function of the magnetic field at the 
physical point. They compare the pressure of a hot pion gas with that of an ideal quark-gluon 
plasma with with a subtracted vacuum energy term due to a nonzero gluon condensate. For weak 
magnetic fields, the transition is first order. The line of first-order transitions ends at critical point 
(y/\qB\, T) = (600, 104) MeV. The critical temperature determined this way is a decreasing function 
of the magnetic field B. 

D'Elia et al have recently carried out lattice simulations in a constant magnetic background [25, 26]. 
They explored various constituent quark masses corresponding to a pion mass of 200 — 480 MeV 
and different magnetic fields, up to \qB\ ~ 20 M 2 . for the lightest quark masses. For these values of 
the pion mass, they found that there is a slight increase in the critical temperature T c for the chiral 
transition. These results have been confirmed by Bali et al [27, 62]. The same group has alo carried 
out lattice simulations for physical values of the pion mass, i.e. M n = 140 MeV. Their results which 
are extrapolated to the continuum limit show that the critical temperature is a decreasing function 
of the magnetic field [27, 62]. Hence the critical temperature for fixed \qB\ as a function of the quark 
mass is nontrivial. This is in stark contrast to most model calculations that imply an increasing 
critical temperature as a functions of B. These include mean-field calculations as well as functional 
renormalization group calculations that incorporate quantum as well as thermal fluctuations [49, 52]. 
Since this discrepancy is not understood, more work is needed to resolve this problem. 
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A Sum-integrals 

In the imaginary-time formalism for thermal field theory, a boson has Euclidean 4-momentum P = 
(Po, p) with P 2 = Pq + p 2 . The Euclidean energy Pq has discrete values: Po = 2n7rT, where n is an 
integer. Loop diagrams involve sums over Po an d integrals over p. With dimensional regularization, 
the integral is generalized to d = 3 — 2e spatial dimensions. We define the dimensionally regularized 
sum-integral by 



and A is an arbitrary momentum scale. The factor (e" /E /4ir) e is introduced so that, after minimal 
subtraction of the poles in e due to ultraviolet divergences, A coincides with the renormalization 
scale of the MS renormalization scheme. 




(A.l) 



where the the integral is defined 




(A.2) 



In the case of particles with electric charge q in a constant magnetic field B, the sum-integral is 
replaced by a sum over Matsubara frequencies Pq = 2irnT, a sum over Landau levels m, and an 
integral over momenta in d — 2 = 1 — 2e dimensions. We then define 
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where the the integral is defined 
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and where the prefactor takes into account the degeneracy of the Landau levels. 
The specific sum-integrals that we need are 
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The functions J n {x) and J^{x) satisfy the recursion relations 
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B Explicit calculations 

In this appendix, we explitcitly calculate the sum-integrals that we need. 



The first sum-integral we need is given by Eq. (A. 5). Summing over Matsubara frequencies we 
can write 
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The first term in Eq. (B.l) is ultraviolet divergent. Calculating it with dimensional regularization, 
we obtain 
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The second term in Eq. (B.l) can be expressed in terms of Jo defined in Eq. (A. 9) using integration 
by parts. Expanding Eq. (B.2) in powers of e through order e, we obtain Eq. (A. 5). We next consider 
the sum-integral (A. 6). Summing over the Matsubara frequencies, we obtain 
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The first integral is ultraviolet divergent and we compute it in dimensional regularization using 
Eq. (B.2) with d = 1 - 2e. This yields 
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Eq. (B.4) shows that the sum over Landau levels m involves the term M B ~ 2e . This sum is diver- 
gent for e = and we regulate it using zeta-function regularization. After scaling out a factor of 
(2\qB\) 1 ~ e , this sum can be written as 
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where x = w^b\ anc ^ £( < ?' s ) = Ylm=o(l + m ) s is the Hurwitz zeta-function. We then find 
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The second term in Eq. (B.3) can be expressed in terms of Jq defined in Eq. (A. 10) using integration 
by parts. Expanding Eq. (B.6) in powers of e, we obtain Eq. (A. 6). 

Taking the derivative of Eq. (A. 5) with respect to M 2 and using the recursion relation (A. 11) for 
J n , one obtains Eq. (A. 7). 

We next consider Eq. (A. 8). It is given by the derivative of Eq. (A. 6) with respect to M B . The 
temperature-independent part of the sum-integral is given by the derivative of Eq. (B.6) with respect 
to M 2 . This yields 

1 2\qB\ fe< E k 2 ' 



2vr ^ 
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Using the recursion relation (A. 12), the temperature-dependent term of Eq. (A. 8) is expressed in 
terms of Jf. 
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